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Abstract 

The discovery of neutrino masses suggests the hkely existence of gauge singlet 
fermions that participate in the neutrino mass generation via the seesaw mech- 
anism. The masses of the corresponding degrees of freedom can range from well 
below the electroweak scale to the Planck scale. If some of the singlet fermions are 
light, the sterile neutrinos appear in the low-energy effective theory. They can play 
an important role in astrophysics and cosmology. In particular, sterile neutrinos 
with masses of several keV can account for cosmological dark matter, which can be 
relatively warm or cold, depending on the production mechanism. The same parti- 
cles can explain the observed velocities of pulsars because of the anisotropy in their 
emission from a cooling neutron star born in a supernova explosion. Decays of the 
relic sterile neutrinos can produce a flux of X-rays that can affect the formation of 
the first stars. Existing and future X-ray telescopes can be used to search for the 
relic sterile neutrinos. 
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1 Introduction 



A discovery in particle physics usually amounts to either a measurement of some parameter 
related to a known particle, or a detection of some new degrees of freedom, new particles. 
One could argue that the discovery of the neutrino mass [1] is likely to be both. Not only is 
it a measurement of the non-zero mass, but it also implies the likely existence of some ad- 
ditional, SU(3)xSU(2)xU(l) singlet fermions, "right-handed" neutrinos. The corresponding 
particles can be made very heavy even for small masses of the active neutrinos (the seesaw 
mechanism [2-6]), but they can also be light, in which case they are called sterile neutrinos. 
The above conclusion is not a mathematical theorem, because the Majorana masses of the 
active neutrinos could emerge as SU(2) triplet higher-dimension operators generated at a 
high scale due to string theory or other dynamics not involving the singlet fermions. How- 
ever, the ease with which the seesaw mechanism can generate these masses makes the seesaw 
mechanism very appealing in its simplicity. This mechanism can explain the smallness of neu- 
trino masses even if the Yukawa couplings are large. However, there is no fundamental theory 
of the Yukawa couplings, and, as discussed below, the naturalness arguments can be made 
in favor of either the large or the small Yukawa couplings. In this review we will consider 
the astrophysical and cosmological ramifications of the seesaw models in which the Yukawa 
couplings and the right-handed Majorana mass terms are allowed to be relatively small. In 
this case, the new degrees of freedom, introduced with the gauge singlet fermions, appear in 
the low-energy effective theory, below the electroweak scale and can have a number of im- 
portant observable consequences, from dark matter to pulsar kicks and other astrophysical 
phenomena. 

As usual, we will call the neutrino weak eigenstates lefthanded if they transform as doublets 
under the SU(2) gauge group of the Standard Model. The gauge singlet fermions in the 
weak eigenstates basis are called singlet or righthanded neutrinos. The mass eigenstates 
corresponding to masses below the electroweak scale will generally have small mixing angles 
and split into two classes: (i) those that couple to Z boson with the coupling approximately 
equal to the electron's coupling to Z, and (ii) those that will have a much smaller coupling 
to Z (suppressed by the mixing parameter smaller than 10~^). The former are called active 
neutrinos^ and the latter are called sterile neutrinos ^ . 



2 Neutrino masses and the emergence of sterile neutrinos 

The Standard Model was originally formulated with massless neutrinos transforming 
as components of the electroweak SU(2) doublets [ot = e,ii,T). To accommodate the 



^ The name sterile neutrino was coined by Bnino Pontccorvo in a seminal paper [7], in which 
he also considered vacuum neutrino oscillations in the laboratory and in astrophysics, the lepton 
number violation, the neutrinoless double beta decay, some rare processes, such as ^ 67, and 
several other questions that have dominated the neutrino physics for the next four decades. 
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neutrino masses, one can add several electroweak singlets Na {a — 1, ...,n) to build a seesaw 
Lagrangian [2-6]: 

C = CsM + iNa^Na - VaaH^ L^N^ - ^ 7V„^7V„ + h.C. (1) 

Here £sm is the Standard Model Lagrangian (with only the left-handed neutrinos and without 
the neutrino masses). We will assume that SU(3)-triplet Higgs bosons [8] are not involved, 
and all the neutrino masses arise from the Lagrangian (1). 



2.1 The seesaw mechanism 



The neutrino mass eigenstates z/}™'' (i = l,...,n + 3) are linear combinations of the weak 
eigenstates {v^, Na}. They are obtained by diagonalizing the (n + 3) x (n + 3) mass matrix: 



^ y^a{H) ^ 

(H) diag{Mi,...,Mjy 



(2) 



As long as all yaa{H) ~ y{H) <^ Ma ~ M, the eigenvalues of this matrix split into two 
groups: the lighter states with masses 

"^(^^1,2:3) ~ ^ (3) 

and the heavier eigenstates with masses of the order of M: 

m(i/('")) - M (a > 3). (4) 

We call the former active neutrinos and the latter sterile neutrinos. Generically, the mixing 
angles in this case are of the order of 

^aa (5) 

but some additional symmetries or accidental cancellations can make them different from 
these generic values. 

One can consider a broad range of values for the number n of sterile neutrinos. Unlike 
the other fermions, the singlets are not subject to any constraint based on the anomaly 
cancellation because these fermions do not couple to the gauge fields. To explain the neutrino 
masses inferred from the atmospheric and solar neutrino experiments, n — 2 singlets are 
sufficient [9], but a greater number is required if the Lagrangian (1) is to explain the r- 
process nucleosynthesis [10], the pulsar kicks [11-18] and the strength of the supernova 
explosion [19,20], as well as dark matter [21-30]. A model often referred to as i/MSM, for 
Minimal Standard Model (MSM) with neutrino (u) masses is the above model with n — 3 
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sterile neutrinos, all of which have masses below the electroweak scale: one has a mass of 
the order of a few kcV, while the two remaining sterile neutrino are assumed to be closely 
degenerate at about 1-10 GeV scale [27]. This model is singled out for the minimal particle 
content consistent with baryogenesis [31,32] and having a dark- matter candidate. However, 
as discussed below, the need for a cosmological mechanism capable of producing colder 
dark matter than that generated by neutrino oscillations may require one to go beyond the 
minimal model, and introduce some additional physics at the electroweak scale [17,33]. 

The scale of the right-handed Major ana masses is unknown; it can be much greater than 
the electroweak scale [2-6] , or it may be as low as a few eV [34] . Theoretical arguments have 
been put forth for various ranges of these Major ana masses. 

2.2 What is more natural? 

The seesaw mechanism [2-6] can explain the smallness of the neutrino masses in the presence 
of the Yukawa couphngs of order one if the Majorana masses are much larger than the 
electroweak scale. Indeed, in this case the masses of the lightest neutrinos are suppressed by 
the ratios {H) /Ma- 

However, the origin of the Yukawa couplings remains unknown, and there is no experimental 
evidence to suggest that these couplings must be of order 1. In fact, the Yukawa couphngs 

of the charged leptons are much smaller than 1. For example, the Yukawa coupling of the 
electron is as small as 10^®. One can ask whether some theoretical models are more likely to 
produce the numbers of order one or much smaller than one. The two possibilities are, in fact, 
realized in two types of theoretical models. If the Yukawa couplings arise as some topological 
intersection numbers in string theory, they are generally expected to be of order one [35], 
although very small couplings are also possible [36]. If the Yukawa couplings arise from the 
overlap of the wavefunctions of fermions located on different branes or in the bulk in extra 
dimensions, they can be exponentially suppressed and are expected to be very small [37-39] . 

In the absence of the fundamental theory, one may hope to gain some insight about the size of 
the Yukawa couplings using 't Hooft's naturalness criterion [40], which states that a number 
can be naturally small if setting it to zero increases the symmetry of the Lagrangian. A small 
breaking of the symmetry is then associated with the small non-zero value of the parameter. 
This naturalness criterion has been applied to a variety of theories; it is, for example, one 
of the main arguments in favor of supersymmetry. (Setting the Higgs mass to zero does not 
increase the symmetry of the Standard Model. Supersymmetry relates the Higgs mass to the 
Higgsino mass, which is protected by the chiral symmetry. Therefore, the light Higgs boson, 
which is not natural in the Standard Model, becomes natural in theories with softly broken 
supersymmetry.) In view of 't Hooft's criterion, the small Majorana mass is natural because 
setting Ma to zero increases the symmetry of the Lagrangian (1) [41,34]. 

One can also ask whether cosmology can provide any clues as to whether the mass scale of 
sterile neutrinos should be above or below the electroweak scale. It is desirable to have a 
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theory that could generate the matter-antimatter asymmetry of the universe. In both hmits, 
of large and small Ma, one can have a successful leptogenesis: in the case of the high-scale 
seesaw, the baryon asymmetry can be generated from the out-of-equilibrium decays of heavy 
neutrinos [42] , while in the case of the low-energy seesaw, the matter- antimatter asymmetry 
can be produced by the neutrino oscillations [31,32]. Big-bang nucleosynthesis (BBN) can 
provide a constraint on the number of light relativistic species in equilibrium [43-46,10,47], 
but the sterile neutrinos with the small mixing angles may never be in equilibrium in the early 
universe, even at the highest temperatures [21]. Indeed, the effective mixing angles of neu- 
trinos at high temperature are suppressed due to the interactions with plasma [48,43,44,49], 
and, therefore, the sterile neutrinos may never thermalize. High-precision measurements of 
the primordial abundances may probe the existence of sterile neutrinos and the lepton asym- 
metry of the universe in the future [50] . 

While many seesaw models assume that the sterile neutrinos have very large masses, which 
makes them unobservable (except, via leptogenesis), it is, of course, worthwhile to consider 
light sterile neutrinos in view of the above arguments, and also because they can explain 
several experimental results, such as dark matter [21], the velocities of pulsars [11,13,15], 
etc. [51-54]. Sterile neutrinos with mass ~ 10 keV and a small mixing angle are predicted 
by some models of neutrino masses [55-57]. In particular, electroweak-singlet fermions with 
masses below the electroweak scale are predicted in some models of dynamical symmetry 
breaking [58-60]. 



2.3 Split scales 



The Majorana masses may, of course, differ in scale. For example, all but one of them 
may be at some high scale, while one Majorana mass may be below the electroweak scale. 
In this case it is useful to integrate out the high-scale states from the full mass matrix in 
eq. (2). Then in the low-energy effective theory one recovers the mass matrix of three active 
neutrinos and one sterile neutrino (see, e.g., Ref. [61] and references therein): 



_^(3X3) _^(1X3)^ 



(6) 



where 



y^(3x3) 



er 



^jue ^jUjU ^/xr 



\ 



Y TTl^e '"-T// '"-TT y 

is the mass matrix of active neutrinos generated by the see-saw mechanism after the high- 
scale states have been integrated out, and Ai^i^^^ = {^ai}, (« = e, yU,r). The patterns of 
masses and mixing angles among the active neutrinos can be affected by the terms induced 
by the presence of sterile neutrinos in the mass matrix of eq. (6). Smirnov and Funchal [61] 



(7) 
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investigated the possible effects of such induced mass terms as a function of sterile neutrino 
mass and the active-sterile mixing. In some ranges of parameters such induced contribution 
could be responsible the observed peculiar features of the effective mass matrix of active 
neutrinos [61]. 

The cosmological bounds on the sum of the active neutrino masses [62] imply that the sum 
of the eigenvalues of mass matrix (6) corresponding to the three (mostly) active eigenstates 
have to be less than 0.54 eV. These eigenvalues, in the limit of small active-sterile mixing, 
are close to the eigenvalues of the matrix (7). Unless some cancellations occur (which is 
possible, especially, in the presence of some new symmetries), one can consider nia < 1 eV. 
One can also assume that Mi ^ m^i and Mi ^ rriap. Then all the active-sterile mixing 
angles are small, and the sterile neutrino has mass close to Mi. In most of our discussion, 
we will concentrate on this one sterile neutrino with mass 

rus^ Ml. (8) 

At the same time, one should keep in mind that some additional light sterile neutrinos may 
exist; they can play an important role in baryogenesis [31,32] and can alter the production 
rates of the dark-matter sterile neutrino in some cases [63]. 

2.4 Can the Majorana masses originate at the electroweak- scale? 

Although the Majorana mass at the keV scale could be a fundamental parameter, it can also 
arise, as the other fermion masses, from the Higgs mechanism [64]. The interaction of the 
form 

Cz^^SNlN^ + h.c, (9) 

where S can be either real or complex, can generate the Majorana mass after the boson 
S acquires a VEV. If S is complex, a U(l) symmetry (lepton number) is spontaneously 
broken by the VEV of S, and the Nambu-Goldstone boson, called Majoron, appears in the 
low-energy effective Lagrangian [64] . This theory does not contain a sterile neutrino that is 
stable enough to be dark matter because the sterile neutrino can decay into the Majoron 
and an active neutrino. If the SU(2)-singlet S is real, then the coupling in eq. (9) breaks the 
lepton number explicitly (just like the mass term in eq. (1)), and no light scalar appears in the 
theory. Hence, sterile neutrino does not have a fast decay channel involving a Majoron, and 
this sterile neutrino is stable enough to be dark matter. The VEV and mass of this scalar will 
be discussed below. A hght, sub-GeV scalar could play the role of the infiaton [65]. However, 
if one requires that the mass and the VEV be of the same order of magnitude, and that 
the dark matter density be acceptable, then one finds that the singlet ends up right at the 
electroweak scale [17,33], where it can have important implications for LHC phenomenology, 
the electroweak phase transition, and baryogenesis [44,66-74]. 

One can consider the following scalar potential: 

V{H, S) = ml\H\^ + mlS^ + X^S^ + Xhs\H\^S^ + X^S^ + X„\H\\ (10) 
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which allows both H and S to have non-zero VEV. After the electroweak symmetry breaking, 
the Higgs doublet and the singlet, each develop a VEV, 



(//) = = 247 GeV, (-5) = -ui. 



(11) 



and the singlet fermions acquire the Majorana masses 



Ma = haVx. 



(12) 



As discussed below, this model is suitable for generating dark matter in the form of sterile 
neutrinos, and it has two very appealing features. First, dark matter generated at the elec- 
troweak scale is not as warm as that produced in oscillations, hence it probably is in better 
agreement with the data. Second, the dark matter abundance, ^dm ~ 0.2 is a natural con- 
sequence of setting the scale of the Higgs singlet, {S) — 10^ GeV, that is at the electroweak 
scale. 



3 Experimental status 

The current limits on sterile neutrino masses and mixing angles are summarized in Figs. 1- 
3. The experimental strategy for discovering the sterile neutrinos depends on the mass and 
mixing parameters. Many experiments have been proposed and performed over the years. Ac- 
celerator and other laboratory experiments are able to set limits or discover sterile neutrinos 
with a large enough mixing angle. The light sterile neutrinos, with masses below 10^ eV and 
large mixing angles, can be discovered in one of the neutrino oscillations experiments [61]. 

In the eV to MeV mass range, the "kinks" in the spectra of beta-decay electrons can be 
used to set limits on sterile neutrinos mixed with the electron neutrinos [75]. Neutrinoless 
double beta decays can probe the Majorana neutrino masses [76]. An interesting proposal is 
to search for sterile neutrinos in beta decays using a complete kinematic reconstruction of 
the final state [77,78]. 

For masses in the MeV-GcV range, peak searches in production of neutrinos provide the 
limits. The massive neutrinos z/j, if they exist, can be produced in meson decays, e.g. tt"^ 
li^Pi, with probabilities that depend on the mixing in the charged current. The energy 
spectrum of muons in such decays should contain monochromatic lines [75] at — (m^ -|- 
m^ — 2m7rm^ — m^.)/2m7r. Also, for the MeV-GeV masses one can set a number of constraints 
based on the decays of the heavy neutrinos into the "visible" particles, which would be 
observable by various detectors. These limits are discussed in Ref. [79]. 

In general, the laboratory experiments, have a limited reach for small mixing angles. As 
discussed below, the flavor mixing angles of sterile neutrinos that can be responsible for the 
pulsar kicks and dark matter are 9 < 10~^. The interaction cross sections of such neutrinos 
are suppresses by factor 6*^ < 10"^^°. The same factor suppresses the oscillation probability. 
Any experiment designed to probe the existence of sterile neutrinos with a small mixing 
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Fig. 1. Experimental and observational limits on sterile neutrinos that have a non-zero mixing 
with the electron neutrino only. This figure is based on the limits from Refs. [79,61,80]. The X-ray 
limits and Lyman-a limits shown here are based on the abundances of relic neutrinos produced by 
neutrino oscillations for zero lepton asymmetry; see discussion in the text and in Refs. [17,81-83] 



angle (see, e.g., Refs. [77,78]) would have to overcome this enormous suppression factor. The 
controversial result from LSND [84-87], later refuted by MiniBooNE [88] had affected much 
of the discussion regarding the sterile neutrinos in the literature. However, the low-mass, 
large mixing angle region probed by these experiments was not motivated at all by particle 
physics, astrophysics, or cosmology [89]. The MiniBooNE results have no bearing on the 
existence of sterile neutrinos in the mass and mixing ranges in which they can account for 
the pulsar kicks and dark matter. 

For masses and mixing angles as small as one needs to explain the pulsar kicks and dark 
matter, the X-ray searches described below probably provide the best opportunity to make 
a discovery. 

Figs. 1-3 show an update of exclusion limits discussed in Refs. [61,79]. We refer the reader to 
these papers for a more detailed discussion and a complete list of references. Here we briefly 
summarize the origin or various limits shown in the figures. 

For all cosmological limits we assume the minimal possible abundance of sterile neutrinos 
consistent with standard cosmology, in which the universe was heated to temperatures above 
1 GeV. (See Ref. [90] for the possibility of very low reheat temperature that could invalidate 
some of the bounds.) Some limits on mass and mixing reported in the literature assume that 
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Fig. 2. Same as in Fig. 1, for the sterile neutrino mixing with only. 

sterile neutrinos constitute 100% of cosmological dark matter, even though their abundance 
would be much lower if no additional production mechanism is introduced besides neutrino 
oscillations. We assume the abundance of sterile neutrinos in halos to be at the lowest value 
consistent with neutrino oscilations. If some additional production mechanisms are at work 
(see, e.g., Refs. [65,17,33]), then the X-ray bounds can exclude some additional parameter 
space as compared to what is shown here. 

The structure formation is affected if sterile neutrinos make a significant contribution to the 
mass density of the universe, i.e., if they make up a non- negligible part of dark matter. The 
bounds marked Large Scale Structure (LSS) are based on the analyzes of Refs. [89,61]. The 
Lyman-a constraints shown in Figs. 1-3 are based on Refs. [81-83]. Part of the reason for the 
difference with Refs. [91-94] is that the authors of Refs. [91-94] assume that sterile neutrinos 
produced by some unspecified mechanism account for the entire observed dark matter, while 
the authors of Ref. [81-83] consider the effects of the amounts of sterile neutrinos that can 
be produced by neutrino oscillations. The oscillations produce a calculable amount of sterile 
neutrinos in standard cosmology, while some additional mechanisms can produce and ad hoc 
amount with a model-dependent free-streaming properties. The latter is a possibility, but it 
is not a good assumption for setting the bounds. 

Likewise, the X-ray limits, based on the non-observation of an X-ray line from the radiative 
decay of sterile neutrinos (see section 7 for discussion) are shown as in Refs. [17,81,95,96], 
based on the actual production rates and not on the assumption that all dark matter is 
accounted for by sterile neutrinos. This makes the excluded region different from that shown 
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Fig. 3. Same as in Fig. 1, for the sterile neutrino mixing with Vr only, 
in Refs. [24,97-113] 

The Cosmic Microwave Background (CMB) bounds are shown as in Refs. [114,61]. The Big 
Bang Nucleosynthesis (BBN) constraints are based on Refs. [115,116,61]. The supernova 
1987A (SN1987A) bounds are based on Refs. [117-121,116,61]. The laboratory bounds, in- 
cluding the accelerator, beta decay, and neutrino oscillations experiments are derived and/or 
summarized in Refs. [79,61], where the reader can find the detailed discussion and references. 



4 Dark Matter in the Form of Sterile Neutrinos 



The singlet fermions are introduced to explain the observed neutrino masses, but the new 
particles can make up the dark matter. Because of the small Yukawa couplings, the keV 
sterile neutrinos are out of equilibrium at high temperatures. They are not produced in the 
freeze-out from equilibrium. However, there are several ways in which the relic population 
of sterile neutrinos could have been produced. 

• Sterile neutrinos could be produced from neutrino oscillations, as proposed by Dodelson 
and Widrow (DW) [21]. If the lepton asymmetry is negligible, this scenario appears to be in 
conflict with a combination of the X-ray bounds [24,97-110,103,111,112,95,96,122,113] and 
the Lyman-a bounds [91,92,123,94]. On the other hand, observations of dwarf spheroids 
point to a non-negligible free-streaming length for dark matter [124-147], which favors 
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warm dark matter, apparently in contradiction with the Lyman-a bounds. It is also pos- 
sible that the sterile neutrinos make up only a fraction of dark matter [81-83], in which 
case they can still be responsible for the observed velocities of pulsars [11,17]. 

• A modification of the DW scenario proposed by Shi and Puller [22] uses a non-zero 
lepton asymmetry L. The oscillations on Mikheev-Smirnov-Wolfenstein (MSW) reso- 
nance [148,149] generate a greater abundance of relic sterile neutrinos with a lower average 
momentum than in the DW case. This results in a colder dark matter with smaller mixing 
angles, which relaxes the bounds from small-scale structure and from the X-ray observa- 
tions. The Shi-Puller (SF) scenario works for a pre-existing lepton asymmetry L > 10~^. 
An economical model that can generate the requisite lepton asymmetry was proposed by 
Laine and Shaposhnikov [63]: decays of the heavier sterile neutrinos could be responsi- 
ble for generating the lepton asymmetry of the universe that creates the conditions for 
producing dark matter in the form of the lighter sterile neutrinos. 

• The bulk of sterile neutrinos could be produced from decays of gauge-singlet Higgs bosons 
at temperatures above the S boson mass, T ~ 100 GeV[17]. In this case, the Lyman-a 
bounds on the sterile neutrino mass are considerably weaker than in DW case or SP case 
because the momenta of the sterile neutrinos are red-shifted as the universe cools down 
fromT ~ 100 GeV [17,33]. 

• Sterile neutrinos could be produced from their coupling to the inflaton [65] or the ra- 
dion [150] . Depending on the time of production, the population of dark- matter particles 
can be warm or cold. Por example if the mass of the inflaton is below 1 GeV [65], one 
docs not expect a significant redshifting of dark matter, which remains warm in this case. 
However, if the sterile neutrinos are produced at a higher scale, they can be red-shifted as 
in the case of the clcctroweak-scale production [17,33]. 

It is important to note that only in the first case, namely, the DW scenario, the dark matter 
abundance is directly related to the mixing angle. In contrast, if the relic population of sterile 
neutrinos arises from the Higgs decays, their abundance is determined by the coupling h in 
eq. (9), while the mixing angle is controlled by a different coupling, as discussed below. We 
also note that both models allow for some production of sterile neutrinos from oscillations, 
but in the case of the singlet Higgs decays (18), the bulk of the sterile dark matter could 
be produced at T ~ 100 GeV, regardless of the value of the mixing angle, which can be 
vanishingly small. Nevertheless, the production by oscillations cannot be turned off, and the 
X-ray bounds, which depend on the mixing angle, apply even in the case when only a fraction 
of dark matter comes from neutrino oscillations [17,81-83]. Of course, the number densities 
of sterile neutrinos in dark halos must be calculated based on the DW production rates. 

Let us now discuss each of these possibilities. 



4-1 Relic sterile neutrinos produced by neutrino oscillations (Dodelson-Widrow scenario) 

Production of sterile neutrinos by oscillations takes place at relatively low temperatures, 
below 1 GeV, and depends only on the mass of the sterile neutrino and its mixing with the 
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active neutrinos. It does not depend the high-scale physics that could be responsible for the 
neutrino masses. Other production mechanisms could be active at some higher temperatures 
and could produce an additional (usually, colder) component of dark matter, as discussed 
below. However, there is no way to turn off the production by oscillations, unless the reheat 
temperature is very low [151,90]. Any additional high-scale mechanism would lead to a two- 
component dark matter. 

Dodelson and Widrow [21] proposed that sterile neutrinos produced by oscillations could 
account for all the dark matter. The calculation of sterile neutrino abundance was subse- 
quently refined to include various corrections and to properly account for the change in the 
QCD degrees of freedom that takes place at around the same time [23,24,26,97,30,152]. 

In general, the singlet neutrino is not an eigenstate of the mass matrix. One can assume, 
for example, that the singlet fermion has a non-zero mixing with the electron neutrino, but 
that the other mixing angles are zero or very small. Then one finds that the mass eigenstates 
have a simple expression in terms of the weak eigenstates: 



|z/|™^) = cos 9m lue) " siu 6*^ \Ni) (13) 
^) = sin^„ \ue) + cosOm \Ni), (14) 

where 9m is the effective mixing angle, which, at finite temperature and density, could be 
different from the vacuum mixing angle. 

If the mixing angle 9m is small, one of the mass eigenstates, z/i behaves very much like a 
pure z/g, while the other, 1^2, is practically "sterile", i.e., its weak interactions are suppressed 
by a factor (sin^ 9m) in the cross section. 

Sterile neutrinos arc produced through oscillations of active neutrinos. The relation between 
their mass and the abundance is very different from what one usually obtains from freeze-out. 
One can trace the production of sterile neutrinos in plasma by solving the corresponding ki- 
netic equations, including the effects of QCD degrees or freedom as discussed in Rcfs. [30,152] 
Here we will give a very crude but intuitive picture of this production mechanism. 

Neutrino oscillations can convert some of the active neutrinos (in equilibrium) into sterile 
neutrinos (out of equilibrium). Matter and the quantum damping effects inhibit neutrino 
oscillations [48,43,44,49]. The mixing of sterile neutrinos with one of the active species in 
plasma can be represented by an effective, density and temperature dependent mixing an- 
gle [21,26,23]: 

3i,2 2^ ^ (AmV2p)^sin^2g 

~ ( Am2 /2p)2 sin^ 29 + ( Am^ /2p cos29 - Vm - V^)^' ^ ' 

Here Vm and Vr are the effective matter and temperature potentials. In the limit of small 
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angles and small lepton asymmetry, the mixing angle can be approximated as 



sin 9. 



(16) 



m 



l+0-27C(l«ifiw)'(g) 



where C = 1-0 for mixing with the electron neutrino, and C = 0.30 for and v^- 

Obviously, thermal effects suppress the effective mixing significantly for temperatures T > 
150 (m/keV)^/^ MeV. If the singlet neutrinos interact only through mixing, all the interaction 
rates are suppressed by the square of the mixing angle, sin^ 6jn- It is easy to see that these 
sterile neutrinos are never in thermal equilibrium in the early universe. Thus, in contrast 
with the case of the active neutrinos, the relic population of sterile neutrinos is not a result 
of a freeze-out ^ . 

In the relevant range of parameters, one can roughly approximate the numerical results for 
the amount of dark matter produced in this scenario [21,23,24,26,97,152]: 



The range of the masses and mixing angles consistent with dark matter and with the X-ray 
bounds discussed below forces the mass of the sterile neutrino to be as low as 1-3 keV. The 
much improved state-of-the-art calculations [30,152] reinforce this conclusion. However, the 
Lyman-a bounds [91,92,123,101,94] appear to disfavor this mass range for the production via 
neutrino oscillations^ . 



4-2 Lepton asymmetry and the Shi-Fuller scenario 

The production scenario proposed by Dodelson and Widrow [21] is altered in the presence 
of a lepton asymmetry of the universe, in which case the production of relic sterile neutri- 
nos can be enhanced by MSW effect [148,149]. Shi and Fuller [22] showed that the MSW 
resonance makes the production more efficient for small missing angles, hence opening up 
some additional parameter space that is less constrained by the X-ray data. In addition, the 
momentum distribution of non-thermal sterile neutrinos produced in this case is colder than 
in the case of zero lepton asymmetry. [22,28,158]. This helps amehorate the tension with the 



^ One immediate consequence of this observation is that the Gershtein-Zeldovich bound [153,154] 
and the Lee-Weinberg bound [155] do not apply to these sterile neutrinos. 

^ Since different production mechanisms can can generate sterile neutrinos with very different free- 
streaming properties for the same mass, the mass bounds from Refs. [92,123,94] do not apply to 
models that consider production by other mechanisms, different from Dodelson- Widrow mecha- 
nism. For example, if sterile dark matter is generated at the electroweak scale, the corresponding 
mass bounds are relaxed by more than factor 3 [17,33,156,157]. 




(17) 
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Lyman-o; bounds ^ . 

For the Shi-Fuller mechanism to work, the Icpton asymmetry should be fairly large, several 
orders of magnitude larger than the baryon asymmetry of the universe. The lepton asym- 
metry of the universe today is not known. Any asymmetry generated at temperatures above 
the electroweak scale would be distributed in comparable amounts between the baryon (B) 
and lepton (L) asymmetries by electroweak sphalerons, which violate (B + L) [159]. How- 
ever, although the standard lore is in favor of much higher temperatures, there is no direct 
observational evidence that the temperature of the universe (after inflation) has at any point 
been higher than the weak scale. In addition, some lepton number violating processes could 
generate a sufficiently large lepton asymmetry at a lower temperature. Affleck-Dine baryoge- 
nesis [160-162] can probably produce very different lepton and baryon asymmetries in some 
cosmological models. 

A very appealing self-contained scenario that allows for the Shi-Fuller production mecha- 
nism was proposed by Laine and Shaposhnikov [63]. They considered the most economical 
model that has enough degrees of freedom to accommodate the active neutrino masses and 
dark matter in the form of sterile neutrinos. This model corresponds to n = 3 in eq. (1) 
and is called i/MSM for Minimal Standard Model with neutrino masses [27]. In addition 
to dark matter, the model is suitable for generating the baryon asymmetry of the universe 
if the two heavier sterile neutrinos have closely degenerate masses M2 ~ M3 ~1-10 GeV 
(with IM2 — M3I ~ 1 keV). The z/MSM leptogenesis [31,32] is similar to the widely studied 
thermal leptogenesis [42] in that a lepton asymmetry L is generated first, by some processes 
involving the heavy singlets, and that this lepton asymmetry is then converted into a baryon 
asymmetry S ~ L. The main difference is that, in thermal leptogenesis, the lepton asym- 
metry is generated by out-of-equilibrium decays of superheavy sterile neutrinos. In the case 
of light sterile neutrinos, the neutrino oscillations replace the decays in generating L ^ 0. 
Oscillations of active neutrinos into sterile neutrinos leave the in-equilibrium plasma with a 
non-zero lepton asymmetry. If the masses of heavy sterile neutrinos, A^2 and in eq. (1) are 
nearly degenerate, the CP violation effects due to some complex phases in the neutrino mass 
matrix arc amplified, thus providing for an acceptable baryon asymmetry of the universe, 
of the order of 10"^'^, consistent with the observations. This scenario is very appealing since 
no additional ingredients are required, besides those that have been added to the standard 
model to accommodate the neutrino masses. ^ 

However, this is not the end of the story. At temperatures well below 100 GeV, when the 
electroweak sphalerons are frozen out, the sterile neutrinos with GeV masses may decay and 
alter the lepton asymmetry of plasma without changing the baryon asymmetry [63] . If this 
happens before the universe reaches the temperatures of the order of 100 MeV, the neutrino 
oscillations responsible for the production of sterile dark matter take place in plasma with a 



See, however, the discussion below regarding some astrophysical evidence in favor of warm dark 
matter. 

^ Although a high degree of mass degeneracy among the heavy sterile states is required, it could 
be a manifestation of some deeper symmetry [163] 
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relatively large lepton asymmetry, which is the pre-requisite for the Shi-Fullcr mechanism. 
Hence, this very economical model can produce both baryon asymmetry and dark matter. 



Although there is still some tension between what's needed for dark matter and the combined 
Lyman-cc bounds and the X-ray bounds, the Shi-Fuller mechanism helps alleviate this tension 
by producing colder dark matter (relevant for small-scale structure) for smaller mixing angles 
(relevant for X-rays). An even colder population of dark matter can be produced by some 
processes at the electroweak scale, as discussed below. 

4-3 Relic sterile neutrinos produced at the electroweak scale 

As discussed above, the Higgs mechanism can give rise to the neutrino Majorana mass if 
the gauge singlet S couples to the neutrinos and acquires a VEV. It is intriguing that, in 
this case, the correct dark matter abundance is obtained for (S) ~ 10^ GeV, suggesting 
that the singlet may be part of an extended Higgs sector [17,33]. Conversely, if one chooses 
the singlet Higgs to have electroweak-scale mass an VEV and if one fixes the coupling as 
required by dark matter, one predicts the sterile neutrino mass in the range where it needs to 
be to explain the pulsar kicks. This production mechanism is worth examining in more detail 
because the sterile neutrinos produced in this scenario make a colder dark matter than those 
produced by oscillations. The reason is entropy production and red shift. The number of 
active degrees of freedom at the electroweak scale is g^{T — 100 GeV) > 110.5, while at the 
low temperature it reduces to (?*(T = 0.1 MeV) = 3.36. The entropy production red-shifts 
the momenta of the dark matter particles and makes this form of dark matter colder, which 
is probably in better agreement with observations. 

Let us consider the following Lagrangian: 



where Co includes the gauge and kinetic terms of the Standard Model, H is the Higgs doublet, 
S is the real boson which is SU(2)-singlet, Lq, (a — e,iJi,T) are the lepton doublets, and Na 
(a = 1, n) are the additional singlet neutrinos. 

The addition of the SNN term opens a new production channel for sterile neutrinos. Dark 
matter particles can be produced in decays S — > A'^A'^ of S particles in equihbrium. At later 
times the sterile neutrinos remain out of equilibrium, while their density and their momenta 
get red-shifted by the expansion of the universe making dark matter "colder" . 

Potential V{H, S) yields the mass of S in the minimum: ms ~ 10^ GeV. One can estimate 
the production of sterile neutrinos by multiplying the S number density (which is ~ for 
T > ms) by the S NN decay rate. 



C^Co + Na {trd^) Na - VaaH L^Na -^S N^Na + V{H, S) + h.C. , 



(18) 



r 



p 



(19) 



(167r) 
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and by the time available for the decay, 

at the latest temperature at which the thermal population of S is still significant, namely 

T ~ ms- (20) 

At lower temperatures, the S number density is too small, much smaller than T^. One obtains 
an approximate result for the number density of dark matter: 



Us 

rp3 



T~ms T"^ 



~ (21) 

r~ms 167r ms 



1 /2 

where Mq = ( ^j^^ ] ~ lO^^GeV is the reduced Planck Mass. This approximate result is 



in good agreement with a more detailed calculation [65,33]. 

Once produced, the dark-matter particles remain out of equilibrium. Production of sterile 
neutrinos via decay S — > NN occurs mainly at temperature ~ ms, after which the 
density of dark matter per co- moving volume undergoes dilution by some factor ^, and the 
particle momenta are redshifted by factor The value of ^ depends on the number of 
active light degrees of freedom at the time of N production. For example, assuming only 
the degrees of freedom in the Lagrangian of eq. (18), that is, the Standard Model with the 
addition of N and S fields, one obtains 

^ «.(r=100GeV) 
^ ^,(r = 0.1MeV) ^ ' 

Based on eq. (21) and taking into account the dilution of the dark matter, it is straightforward 
to show that the S boson decays produce 

--(i^)'O(f)- 



On theoretical grounds, it is reasonable to expect a gauge singlet in the Higgs sector; such 
an extension of the standard model has been the subject of many studies [164,67-70,17,71- 
74,33]. For the electroweak-scale mass and VEV, 

(S) r^ms^ 10^ GeV, 

the Yukawa coupling 

/ ~ io-« 

is required to generate the correct amount of dark matter (see eq. (23)). So far we just fit one 
parameter by choosing the value of another parameter. However, this value of the Yukawa 
coupling turns out to be very special. The corresponding Majorana mass M is 

Ms = f{S) ~ 10-^ X (10^ - 10^) GeV ~ (1 - 10) keV. (24) 
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This mass is exactly what is needed to explain the pulsar velocities by the anisotropic 
emission of sterile neutrinos, as discussed below. If the value of fig in eq.(23) was tuned, the 
mass in eq.(24) is the outcome of the calculation. This is intriguing and encouraging. One 
can see this as an extra motivation for the search in the relevant mass range. 

Not only is (1-20) keV mass good for the pulsar kicks, it is also in the most interesting range 
for dark matter. Indeed, if the mass came out to be lower, it would violate the Tremaine- 
Gunn bound [165] and its generalizations [131,156,157,166,167]. If the mass came out to be 
much greater, the sterile neutrino would decay too fast to be dark matter. Therefore, the 
mass range (1-20) keV is favored by both the pulsar kicks and dark matter. We note that 
theoretical models of neutrino masses can readily produce a sterile neutrino with a required 
mass [168-171]. 

The dilution of dark matter by entropy production affects the average momentum of the dark- 
matter particles. In the case of production at the electroweak scale, the average momentum 
of the dark matter particles, relative to the temperature, is 



(Ps)(T«:lMeV) = 0.76T 



110.5 

,^*(^s). 



1/3 



(25) 



We have derived eq. (21) and eq. (23) using some very simple estimates of dark matter 
production at the electroweak scale. This admittedly simplistic treatment is in good agree- 
ment with the exact solution based on the corresponding kinetic equations [65,17,33]. The 
momentum distribution of the resulting dark matter was studied in detail in Refs. [156,157]. 

The Lagrangian in eq.(18) can describe some very interesting Higgs dynamics, which can 
have a first-order phase transition [67,71,74]. In this case the standard electroweak baryo- 
genesis [159] can take place in the course of this first-order phase transition. The model of 
eq. (18) can easily be modified to include a sufficient amount of CP violation: all that is 
required for a successful baryogenesis is to include the second Higgs doublet [67]. 

The first-order phase transition can, in principle, further redshift the dark matter if the dark 
matter particles are produced mostly before the transition. This possibility was analyzed 
in detail in Ref. [33]. Finally, one should also consider the possibility of S decays out of 
equilibrium, which can occur if the doublet-singlet couplings are small. This possibility was 
also considered in Ref. [33]; it does not result in a significant production of dark-matter 
particles for generic ranges of model parameters. 



4-4 How cold is dark matter in the form of sterile neutrinos? 



While the existence of dark matter is well established by a consensus of several indepen- 
dent lines of evidence [172], the data do not allow to discriminate between a wide variety of 
candidate particles [172,173]. The list of such particle candidates includes the hghtest super- 
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symmetric particles (LSP) and Q-balls predicted by super symmetry, axions, sterile neutrinos, 
and a variety of other candidates, well motivated from the theoretical point of view. It would 
be a major breakthrough, of course, to obtain some additional information about dark mat- 
ter directly from astrophysical data. One possibihty to learn about dark matter is to study 
the small-scale structure, where the particle candidates give very distinct predictions. 

Expressions such as "cold dark matter" (CDM) and "warm dark matter" (WDM) are often 
used, but they do not capture the broad variety of small-scale properties exhibited by different 
dark matter candidates. A long hst of particle candidates [172] labeled traditionally as CDM 
or WDM fit the observed structure equally well on large scales. However, on small scales, 
their predictions differ. The smallest structures in the universe are directly related to particle 
properties and production history of dark matter. For instance, if all dark matter is in the 
form of axions, then the small-scale cutoff in the matter power spectrum is as small as 
10~^^M^ [174]. Supersymmetric dark matter in the form of LSP has a much larger cutoff, 
(10~^ — 1O^)M0 [175], and the dark matter in the form of sterile neutrinos produces an even 
larger cutoff [17,33,156,157], which depends on the mechanism by which these particles were 
produced [17]. For the sake of unambiguous terminology, one can define CDM as idealized 
dark matter that has no small-scale suppression in the matter power spectrum, and we 
define WDM as dark matter that has a non-negligible value of such a cutoff. A definitive 
determination of the small-scale power spectrum would be a powerful discriminator between 
different candidates. 

The data regarding the structures on small scales are becoming available but the observations 
are often difficult to interpret. On the one hand, there are observations of the Lyman-a forest 
from gas clouds at high redshifts, which are interpreted, through modeling, as evidence of 
small-scale structure [92.123,176.177,94,82,83]. On the other hand, there are observations 
of the motions of stars in dwarf spheroid galaxies, which yield the gravitational potentials 
created by dark matter in these dark-matter dominated objects [125,136-139,141-147]. These 
observations are consistent with cored profiles that are most easily accommodated in the case 
of warm dark matter with a suitable small-scale cutoff in the matter power spectrum. 

Claims have been made that CDM suffers from several inconsistencies between the N-body 
simulations and the observations. These inconsistencies may be resolved by better modeling, 
but at present they present a challenge. It is intriguing, however, that WDM is free from 
many of these problems altogether. One challenge to CDM is a discrepancy between the 
number of satellites predicted in N-body simulations and the number of observed satellites 
in galaxies such as the Milky Way [124,127,129,128,178]. Warm dark matter (WDM) can 
suppresses small scale structure and dwarf galaxy formation [130], hence ameliorating this 
discrepancy. Several other potential problems in the CDM paradigm may be solved by the 
reduction of power on small scales due to WDM. For example, WDM can reduce the number 
of halos in low-density voids bringing it in better agreement with observations [132,130]. 
Another problem is the relatively low densities of the galactic cores implied by the rotation 
curves as compared to what is predicted from the ACDM power spectrum [131]. This problem 
can also be rectified by a reduction of the initial power spectrum of density fiuctuations on 
small scales [133]. There appears to be a so called "angular-momentum" problem of CDM 
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halos, in which the gas cools at very early times into small mass halos and leads to massive 
low-angular momentum gas cores in galaxies [126]. Another apparent inconsistency is the 
lack of disk-dominated or pure-disk galaxies predicted in CDM models [179] and the high 
degree of correlations and organization seen in the statistical properties of galaxies [180,181]. 
Both of these inconsistencies can be eliminated if the merger rates were lower than expected 
in the case of CDM. All of these inconsistencies arise in CDM models on small scales. In our 
opinion, none of these challenges can rule out CDM: a better understanding of the N-body 
dynamics, the galaxy formation history and the possible feedback mechanisms may provide a 
resolution for each of these apparent inconsistencies [182] . It is true, however, that warm dark 
matter in the form of sterile neutrinos is free from all these small-scale problems altogether, 
while on large scales WDM fits the data as well as CDM. From the particle physics point of 
view, we see no reason to favor one or the other possibility, because well-motivated candidates 
exist in both cases. Improved understanding of structure formation and assembly history in 
the case of warm dark matter presents a number of interesting observational consequences 
that can be tested by future observations, such as, for example, the filamentary structure of 
halos ^ expected from WDM [183]. 

The momentum distribution of dark matter particles determines how warm the dark matter 
is. This, in turn, depends on the production mechanism [131,17,156,157]. A very useful 
parameter is the phase density that can be defined using the dark matter density p and the 
velocity dispersion {v^): 

Observations of dwarf spheroidal galaxies can be used to infer bounds on the primordial value 
of Q. These data have been interpreted as implying a ~ 1 keV mass for sterile neutrinos, if 
these particles make up all of dark matter [157,184]. 

One can also compare different production mechanisms mechanisms using the free-streaming 
length: 

The average momentum (p^) depends on the way the sterile dark matter is generated in the 
early universe: 



0.9 for production off — resonance 

0.6 for MSW resonance (depending on L) (28) 
0.2 for production at T > 100 GeV 



This is an admittedly simplistic way to present the results for the small-scale power spec- 
trum. However, it allows one to juxtapose the three different scenarios and to illustrate their 
differences. One can see that the three scenarios produce dark matter that becomes "warm" 

^ This could affect the Lyman-a bounds. 
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on different length scales. Hence, one can hope that the future astronomical observations 
will be able to distinguish between such scenarios. 

The clustering properties of dark matter in the form of sterile neutrinos have been studied in 
detail [97,185-188,140,156,189,157,166,167], and these studies conclude, in particular, that, if 
the population of dark-matter sterile neutrinos originate at or above the electroweak scale, as 
in the model of Refs. [17,33], then their power spectrum is consistent with the data on small- 
scale structure, while the other parameters, such as mass and the mixing angle, are within 
their allowed ranges [156,189,157]. The phase space arguments [165,131,156,157,166,167] can 
be used set the lower bound on the mass of dark-matter sterile neutrinos, which must be 
heavier than 0.5 keV if they make up all the dark matter and if they were produced at the 
electroweak scale, with a subsequent cooling, as in Refs. [17,33,156,157]. 



5 Sterile neutrinos and the supernova explosion 

Sterile neutrinos with masses below a GeV can be produced in supernova explosions. They 
can play an important role in the nucleosynthesis [10,190], as well as in generating the 
supernova asymmetries and the pulsar kicks [11]. The role played by the sterile neutrinos in 
the supernova depends on the particle mass. 

5.1 Pulsar velocities 

Radio pulsars are magnetized neutron stars. They have some very high velocities, the origin 
of which remains unknown. The velocities of pulsars are measured either by observation 
of their angular proper motions [191-194], or by measuring the velocity of an interstellar 
scintillation pattern as it sweeps across the Earth [195-199] . Each of the two methods has 
certain advantages. Using the former method, one can get very precise measurements with 
the help of a high-resolution radio interferometer, but such observations take a long time. 
Measuring the velocity of a scintillation pattern can be done quickly, but the inference of 
the actual pulsar velocity must rely on some assumptions about he distribution of scattering 
material along the line of sight. (For instance, if the density of scatterers is higher near 
Earth, the pattern moves slower than the pulsar.) In addition to observational errors, one 
has to take into account various selection effects. For example, fast and faint pulsars are 
under-represented in the data as compared with the slow and bright ones. Therefore, one 
has to carefully model the pulsar population to calculate the three-dimensional distribution 
of pulsar velocities corresponding to the observed two-dimensional projection of their proper 
motion [200-203]. 

Based on the data and population models, the average velocity estimates range from 250 km/s 
to 500 km/s [195-202]. The distribution of velocities is non-gaussian, and there is a substan- 
tial population of pulsars with velocities in excess of 700 km/s. Some 15% of pulsars [201,202] 
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appear to have velocities greater than 1000 km/s, while the fastest pulsars have speeds as 
high as 1600 km/s. Obviously, an acceptable mechanism for the pulsar kicks must be able 
to explain these very fast moving pulsars. In addition, one hopes that the pulsar spins can 
be explained by the same kick mechanism [204]. 

Pulsars are born in supernova explosions, so it would be natural to look for an explanation 
in the dynamics of the supernova [205-207]. However, 3-dimensional numerical calculations 
often fail to explain pulsar kicks in excess of 200 km/s [206]. The exciting recent discovery of 
standing accretion shock instabihty (SASI) [208,209] could lead to an explanation of both the 
pulsar kick and the spin. However, it appears that this mechanism predicts a misalignment 
between the axis of rotation and the direction of the pulsar velocity because the spiral SASI 
mode, instrumental in spinning up the pulsar, is expected to give it a kick in the direction 
orthogonal to the axis of rotation [210]. At the same time, the observational evidence is 
growing in favor of alignment of the direction of motion with the axis of rotation [211-219]. 
(As discussed below, such an alignment is a generic prediction of the neutrino driven kick 
mechanisms.) 

The hydrodynamic kick could be stronger if some large initial asymmetries developed in 
the cores of supernova progenitor stars prior to their collapse. Goldreich et al. [220] have 
suggested that unstable g-modes trapped in the iron core by the convective burning layers 
and excited by the e-mechanism may provide the requisite asymmetries. However, according 
to numerical calculations [221], the e-mechanism may not have enough time to significantly 
amplify the g-modes prior to the collapse. Evolution of close binaries [222] and asymmetric 
emission of radio waves [223] have been considered as possible causes of the rapid pulsar 
motions. However, both of these explanations fail to produce a large enough effect. 

Most of the supernova energy, as much as 99% of the total 10^^ erg are emitted in neutrinos. A 
few per cent anisotropy in the distribution of these neutrinos would be sufficient to explain the 
pulsar kicks. This could be an alternative to a hydrodynamic mechanism which would require 
a much larger asymmetry in what remains after the neutrinos are subtracted from the energy 
balance. The numerical calculations of the supernova assume that neutrino distribution is 
isotropic. What if this is not true? 

Since the total energy released in supernova neutrinos is £■ ~ 3 x lO^^erg, the outgoing 
neutrinos carry the total momentum 



Pi.,totai ~ 1 X 10*^^g cm/s. 



(29) 



A neutron star with mass IAMq and speed v — 1000 km/s has momentum 



= (1.4Mo)vRi3 X 10' 
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A few per cent asymmetry in the neutrino distribution is, therefore, sufficient to explain the 
observed pulsar velocities. What could cause such an asymmetry? The obvious suspect is the 
magnetic field, which can break the spherical symmetry and which is known to have an effect 
on weak interactions. We will examine this possibility in detail. Rather than discussing the 
calculations presented in the original papers, we will follow a simple (and more pedagogical) 
discussion from Ref. [15]. 

5.2 Why a sterile neutrino can give the pulsar a kick 

Active neutrinos are always produced with an anisotropy in a strong magnetic field, but they 
usually escape isotropically. The asymmetry in production comes from the asymmetry in the 
basic weak interactions in the presence of a magnetic field. Indeed, if the electrons and other 
fermions are polarized by the magnetic field, the cross sections of the urea processes, such as 

n + e'^ ^ p + Pe: 

p + e~ ^n + Vei 

depend on the orientation of the neutrino momentum with respect to the electron spin. The 
polarization of electrons creates an asymmetry: 

(7(Te-,Ti^)7^(7(Te-,ii/) (31) 

Depending on the fraction of the electrons in the lowest Landau level, this asymmetry can be 
as large as 30%, which is, seemingly, more than one needs to explain the pulsar kicks [224]. 
However, this asymmetry is completely washed out by scattering of neutrinos on their way 
out of the star [225,226]. This is intuitively clear b result of scatterings, the 

neutrino momentum is transferred to and shared by the neutrons. The neutrinos undergo 
multiple scattering and remain almost in equilibrium as they diffuse out of the protoneutron 
star. In the approximate thermal equilibrium, no asymmetry in the production or scattering 
amplitudes can result in a macroscopic momentum anisotropy. This statement can be proved 
rigorously [225,226]. 

However, if the neutron star cooling produced some particles whose interactions with nuclear 
matter were even weaker than those of ordinary neutrinos, such particles could escape the 
star with an anisotropy equal the anisotropy in their production. The state v'^^ in equation 
(14), whose interactions are suppressed by (sin^ ^^) can play such a role. This could be the 
explanation of the observed velocities of pulsars. It is intriguing that the same particle can 
make up the dark mater. 

Let us again consider a model with only one sterile neutrino in which the mass eigenstates are 
admixtures of active and sterile neutrinos, as in equation (14). For a sufficiently small mixing 
angle in matter, ^m, between z/g and z/^, only one of the two mass eigenstates, is trapped 
in the core of a neutron star. The orthogonal state, i'2^\ escapes from the star freely. This 
state is produced in the same basic urea reactions {e.g., n + e"*" ^ P + ^e ^-nd p + e" ^ n-\-Vg) 
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with the effective Lagrangian couphng equal the weak couphng times sin 9m- The production 
and anisotropy can be greatly enhanced if the active neutrinos undergo a resonant conversion 
into the sterile neutrinos at some density [11]. This effect plays an important role in some 
range of parameters, although this kind of enhancement is not a necessary condition for the 
pulsar kick. We will consider two ranges of parameters, for which the Ue fs oscillations 
occur on and off resonance. 

5.3 Pulsar kicks from active-sterile neutrino oscillations 

Let us consider neutrino cooling during the first 10-15 seconds after the formation of a hot 
protoneutron star. Depending on the mass and the mixing angle, there may or may not be a 
resonant conversion of the active to sterile neutrinos at some density in a hot neutron star. 
If there is an MSW resonance, the position of the resonance point depends on the density 
and the magnetic field. The latter introduces the required anisotropy. In the absence of the 
MSW resonance, an off-resonance emission from the entire volume of the neutron star core is 
possible. We will see that this emission is efficient only after the matter potential has evolved 
from its initial value to nearly zero. This important evolution [23] requires some time, and 
can cause a delayed kick [18], which, in turn, imposes constraints on the masses and mixing 
angles. We will consider the following three possibilities for the pulsar kick: 

(1) MSW resonance in the core (p> 10^4 g/cm^) 

(2) MSW resonance outside the core (p< lO^^g/cm^) 

(3) an off-resonance emission from the core 

The three regimes are probably mutually exclusive. For example, for all the masses that 
are consistent with the resonance, the matter potential evolves very slowly, and there is no 
significant emission from the core off-resonance. 

5.4 MSW resonance in the core 

For simplicity, let us assume that inside the star there exists a uniform (dipole) magnetic 

— * 

field B. Neutrino oscillations in a magnetized medium are described by an effective potential 
that depends on the magnetic field [227-235] in the following way: 

yK)=0, (32) 

v{v,) = -y(Pe) = 1^0 (31; - 1 + 4y.j, (33) 

V{v,,r) = -V{u,,r) = (n - 1 + 2 + ^ (^j (34) 

where Fg O^v^) is the ratio of the number density of electrons (neutrinos) to that of neutrons, 
S is the magnetic field, ^ is the neutrino momentum, Vq — 10eV(p/10^^g/cm^). The magnetic 



25 



field dependent term in equation (34) arises from polarization of electrons and not from a 
neutrino magnetic moment, which in the Standard Model is small, and which we will neglect. 
(A large neutrino magnetic moment can produce a pulsar kick through a different mechanism 
proposed by Voloshin [236].) 

The condition for a resonant MSW conversion z/j ^ uj is 



2k 



cos 2% + V{ui) 



2k 



cos 29ij + V{vj) 



(35) 



where fjj- can be either a neutrino or an anti-neutrino. 



In the presence of the magnetic field, condition (35) is satisfied at different distances r± 
from the center (Fig. 4), depending on the value of the {k ■ B) term in (34). The average 
momentum carried away by the neutrinos depends on the temperature of the region from 
which they escape. The deeper inside the star, the higher is the temperature during the 
neutrino cooling phase. Therefore, neutrinos coming out in different directions carry momenta 
which depend on the relative orientation of k and B. This causes the asymmetry in the 
momentum distribution. 
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Fig. 4. For MSW resonance in the core, the sterile neutrino energy depends on the temperature 
around the resonance point. 



The surface of the resonance points is 

r(0) = ro + 5 COS0, (36) 

where COS = {k-B)/{kB) and 5 is determined by the equation {dNn{r)/dr)6 ~ e(3Ai'e/vr^) B. 
This yields 



5=^5/^, (37) 

where /Xg ~ (37r^A^e)^^'^ is the chemical potential of the degenerate (relativistic) electron gas. 

In the core of the neutron star, at densities above 10^^ g/cm^, one can assume the black-body 
radiation luminosity in sterile neutrinos: 

F,^(r)oc T^(r). (38) 



26 



Then the asymmetry in the momentum distribution is 

Ak, ^ 1 T\r + S)-T^{r-6) ^AldT 

~ 3 THr) 3T^^^^^' ^^^^ 

where a factor (1/3) represents the result of integrating over angles. 

Now we use the expression for 5 from eq. (37) and replace the ratio of derivatives (^) / (^^) 
bv — • 

ks ~37r2 \TdNj^- 
To calculate the derivative in (40), we assume approximate thermal equilibrium. Then one 
can use the relation between the density and the temperature of a non-relativistic Fermi gas: 

_ 2(m,T)3/2 |- ^zdz 



in I v^^^ (41) 



%/2- 

where m„ and /x^ are the neutron mass and chemical potential. The derivative (dT/dNn) 
can be computed from (41). Finally, 

A^^W2^^_ (42) 

We have assumed that only one of the neutrino species undergoes a resonance transition into 
a sterile neutrino. The energy, however, is shared between 6 species of active neutrinos and 
antineutrinos. Therefore, the final asymmetry due to anisotropic emission of sterile neutrinos 
is 6 times smaller: 
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k 6 ks 37r2 ^3/2^2 

= 0.01 (^^) fn^y^' C^SlV ( ^ ) (43) 
VlOOMeVy VSOMeVy \ T J \3x10^^gJ ^ ' 

This estimate [11] can be improved by considering a more detailed model for the neutrino 
transport and by taking into account time evolution of chemical potentials discussed below. 
However, it is clear that the magnetic field inside the neutron star should be of the order 
of 10^^ G. The approximation used in equation (38) holds as long as the resonant transition 
occurs deep in the core, at density of order 10^'^gcm~3. This, in turn, means that the sterile 
neutrino mass must be in the keV range. 



5.5 Resonance at densities below 10^'^ g/cm^ 



For smaller masses, the resonance occurs at smaller densities. Outside the core, fewer neu- 
trinos are produced, while there is a flux of neutrinos diffusing out of the core. Therefore, 
the approximation (38) is not valid. 
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Outside the core, the active neutrinos can interact with matter and deposit momentum to 
the neutron star medium. After an active neutrino is converted into a sterile neutrino, it 
no longer interacts with matter and comes out of the star. Some of the electron neutrinos 
are absorbed on their passage through a layer in the neutron star atmosphere. In particular, 
there is a charged-current process 

Vf,n — > e~p^ , (44) 

which has a cross section a = 1.8 G"^ E^, where E,y is the neutrino energy. If the resonant 
conversion z/g — ^ Ug occurs at different depths for different directions, the neutrinos may 
spend more time as active on one side of the star than on the other side of the star, as shown 
in Fig. 5. Hence, they deposit more momentum through their interactions with matter on 
one side than on the other side. Let us estimate this difference. (This argument was used 
before in application to active neutrino transport [12]. Here we adopt it to sterile neutrinos.) 

Depending on the magnetic field, the resonance lies at different depths, eq. (36). Hence, the 
neutrinos on one side of the star pass an extra layer of thickness (25 cos 0) as active, while 
the neutrinos on the other side pass this layer as sterile. The active neutrinos going through 
a layer of nuclear matter with thickness 2S have an extra probability 

Ps=i26)crNn 

to interact and deposit momentum k ^ to the neutron star. This momentum is not 
balanced by the neutrinos on the other side of the star because they go through the this 
layer as sterile neutrinos. The cross section a ~ l.'^G^pE'^. 
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Fig. 5. For MSW resonance outside the core, the neutrino passes between r_ = ro — 6cos(j) and 
r_|_ = ro + 5 cos as a sterile Vs on one side of the star, while it still propagates as an active Va on 
the other side. The active neutrinos interact and deposit some extra momentum on the right-hand 
side, between r_ and r+. Since the neutron star is a gravitationally bound object, the momentum 
deposited asymmetrically in its outer layers gives the whole star a kick. 

Obviously, the neutron star as a whole is a gravitationally bound object, so any momentum 
deposited on one side of the star gives the whole neutron star a kick. 

The difference in the momentum deposition per electron neutrino between the directions 
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and —0 is 

= (25 cos (/)) 7V„ a = 1.8 GlEl ^ ^ hj,^ cos (45) 

where we have used eq. (37) and introduced the scale height of the electron density — 
[dOnNe)/dr]-\ 

We take ^ 0.1, ^ 3T ^ 10 MeV, fi^ ~ 50 MeV, and h^^ ^ 6 km. We assume 
T a; 3 MeV because it is a realistic average temperature around the neutrinosphere, in 
agreement with theoretical models as well as observations of the supernova SN1987A [117]. 
(This temperature is lower than the core temperature used earlier.) 

After integrating over angles and taking into account that only one neutrino species under- 
goes the conversion, we obtain the final result for the asymmetry in the momentum deposited 
by the neutrinos. 



k °'°^(3Mev) (50Mev) (ekm) (iQISg) ' ^^^^ 

This is, clearly, a sufficient asymmetry for the pulsar kick. The corresponding region of 
parameters is shown as region "2" in Fig. 13. 

The above estimates are valid as long as 6 is much smaller than the mean free path. One 
can also describe the propagation of neutrinos in this region using the so called diffusion 
approximation. It was used for the neutrino transport by Schinder and Shapiro [237] in 
planar approximation, and it can be applied to the pulsar kicks [12,238,239,14]. 



5.6 Off-resonance transitions 



Let us now consider the case of the off-resonance emission from the core. For masses of a 
few keV, the resonant condition is not satisfied anywhere in the core. In this case, however, 
the off-resonant production of sterile neutrinos in the core can occur through ordinary urea 
processes. A weak-eigenstate neutrino has a sin^ 9 admixture of a heavy mass eigenstate 
z/2. Hence, these heavy neutrinos can be produced in weak processes with a cross section 
suppressed by sin^ 9. 

Of course, the mixing angle in matter 9„i is not the same as it is in vacuum, and initially 
sin^^^„i ^ sin^ However, as Abazajian et al. [23] have pointed out, in the presence of 
sterile neutrinos the mixing angle in matter quickly evolves toward its vacuum value. When 
sin^ 9m ~ sin^ 9, the production of sterile neutrinos is no longer suppressed, and they can 
take a fraction of energy out of a neutron star. 

Following Abazajian, Fuller, and Patel,[23] one can estimate the time it takes for the matter 
potential to evolve to zero from its initial value V^^^Vg) ~ (—0.2... -|- 0.5)Vo. The time scale 
for this change to occur through neutrino oscillations off-resonance is 
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^off-res., 4v/27r^^. (y(0)(^.))3 1 

6x IQ-'^g / y(°)(i/e) y /SOMeVy /lOkeV^y 

As long as this time is much smaller than 10 seconds, the mixing angle in matter approaches 
its vacuum value in time for the sterile neutrinos to take out some fraction of energy from a 
cooling neutron star. 

The urea processes produce ordinary neutrinos with some asymmetry depending on the 
magnetic field [224]. The same asymmetry is present in the production cross sections of 
sterile neutrinos. However, unlike the active neutrinos, sterile neutrinos escape from the star 
without rescattering. Therefore, the asymmetry in their emission is not washed out as it is 
in the case of the active neutrinos [225,226]. Instead, the asymmetry in emission is equal the 
asymmetry in production. 

The number of neutrinos dN emitted into a solid angle dVl can be written as 

dN 

— = iVo(l + ecose,), (48) 

where is the angle between the direction of the magnetic field and the neutrino momen- 
tum, and A^o is some normalization factor. The asymmetry parameter e is equal 

^ = 4-A*o i^) , (49) 

9i + m \EtotJ 

where and are the axial and vector couplings, E'tot and Eg are the total neutrino 
energy and the energy emitted in sterile neutrinos, respectively. The number of electrons in 
the lowest Landau level, ko, depends on the magnetic field and the chemical potential /i, but 
it can be in the 0.1 < ko < 0.7 range for some reahstic values of these parameters [13]. 

The momentum asymmetry in the neutrino emission is 



where is the fraction of energy carried by the sterile neutrinos. To satisfy the constraint 
based on the observation of neutrinos from supernova SN1987A, we require that < 0.7. 
As can be seen from Fig. 6, the asymmetry in equation (50) can be of the order of a few per 
cent, as required, for magnetic fields 10^^ — 10^^ G. 

Surface magnetic fields of ordinary radio pulsars are estimated to be of the order of 10^^ — 
lO^^G. However, the magnetic field inside a neutron star may be much higher [240-242], 
probably up to lO^^G. The existence of such a strong magnetic field is suggested by the 
dynamics of formation of the neutron stars, as well as by the stability of the poloidal magnetic 
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Fig. 6. The fraction of electrons in the lowest Landau level as a function chemical potential. The 
value of the magnetic field is shown next to each curve [13]. 

field outside the pulsar [241]. Moreover, the discovery of soft gamma repeaters and their 
identification as magnetars [241,242], i.e., neutron stars with surface magnetic fields as large 
as 10^^ G, gives one a strong reason to believe that the interiors of many neutron stars may 
have magnetic fields as large as 10^^ — 10^^ G. There are also plausible physical mechanisms 
that can generate such a large magnetic field inside a cooling neutron star [241]. The a — Q 
dynamo mechanism makes the fields grow during the first 10-15 seconds after the formation 
of the neutron star [241]. The corresponding equations can be analyzed in linear regime, 
until the perturbations grow large. The saturation field of the order of 10^^ G, at which the 
linear analysis is no longer valid is probably a good estimate of the interior field. 

5.7 Pulsar kicks from the active neutrinos alone? 

One can ask whether the sterile neutrino is necessary and whether the oscillations of active 
neutrinos alone could explain the pulsar kicks. The interactions of muon and tau neutrinos 
in nuclear matter are characterized by a smaller cross section than those of the electron 
neutrinos. This is because the electron neutrinos z/g interact through both charged and neutral 
currents with electrons, while z/^ and Ur interact with electrons through neutral currents 
alone. Therefore, nuclear matter is less transparent to z/g than to z/^^r, ^/^,t- As a result, the 
surface of last scattering for z/^,- and z/^^-r lies (about a kilometer) deeper than that of Ug. 
The electron antineutrino can interact through charged currents with positrons while they 
are present in nuclear matter. The z/g mean free path starts out closer to that of z/g, but, as 
the number of positrons diminishes during the cooling period, this mean free path increases 
and becomes comparable to that of z/^,,- and u^^r- 

Since the fi— and r— neutrinospheres lie inside the electron neutrinosphere, it is possible 
that neutrino oscillations could convert a z/g into z/^ or z/^ at some point between the two 
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neutrinospheres, where the Ue is trapped, but the i/^^^- is free-streaming. Then the shift in 
position of the MSW resonance would result in an anisotropy of the outgoing momentum. 
This mechanism could explain the pulsar kicks, but it would require one of the active neutrino 
masses to be of order 10^ eV [243]. This is not consistent with the present experimental data 
on neutrino masses. 



5.8 What if neutrinos have a large magnetic moment? 

The neutrino magnetic moment in the Standard Model is very small, iiy ~ 3xlO^^^(mi//eV)/ij,, 
where ii^ = e/2me and rrii, is the (Dirac) neutrino mass. This is why we have so far neglected 
any effects of direct neutrino interactions with the magnetic field. 

However, the present experimental bounds allow the neutrino magnetic moments to be as 
large as 10~^^iig. If, due to some new physics, the neutrino magnetic moment is large, it may 
open new possibilities for the pulsar kick. Voloshin [236] has proposed an explanation of the 
pulsar kick based on the resonant spin-flip and conversion of the left-handed neutrinos into 
the right-handed neutrinos, which then come out of the neutron star. Voloshin argued that 
the magnetic field inside a neutron star may be irregular and may have some asymmetrically 
distributed "windows" , through which the right-handed neutrinos could escape. The resulting 
asymmetry may, indeed, explain the pulsar velocities. 

Other kinds of new physics may cause the pulsar kicks as well [244,238,245]. 



5.9 Time delays 

As one can see from eq. (47), there can be a non-negligible time delay prior to the onset 
of an efficient emission of sterile neutrinos. The matter potential relaxes to zero on some 
time scales related to the mass and mixing parameters of the sterile neutrinos [18]. After the 
equilibrium is achieved at Vm ~ 0, the effective mixing angle in matter is close to that in 
vacuum, and the emission of sterile neutrinos proceeds at a much higher rate. There is a con- 
siderable uncertainty in the equilibration time given by eq. (47) because several parameters 
are functions of time and position in the star and are not known precisely. In addition, the 
equilibration does not have to occur simultaneously in the entire star. Nevertheless, the very 
possibility of the delay and its connection with the fundamental neutrino parameters can 
allow one to distinguish this mechanism from some alternatives. The studies of the pulsar 
populations can eventually determine the delay based on the observational parameters [218]. 
The allowed ranges of parameters corresponding to some time delays are shown in Fig. 7. 



32 




20.0 



2.0 



delayed pulsar kick 
(no MSW) 



1.5 



1.0 



1x10 



.-II 



2x 10 



5x 10 



1x10 



,-10 



2x 10 



,-10 



5x10 



,-10 



1x10 



sin^e 



Fig. 7. The allowed regions for delayed kicks with delays from 1 through 5 seconds (assuming the 
other parameters are fixed) are shown by black solid lines marked by the numbers representing the 
delay time in seconds [18]. 

5.10 B — V correlation? 

The neutrino kick mechanism does not predict a correlation between the direction of the 
surface magnetic fields and the pulsar velocity. The kick velocity is determined by the mag- 
netic field inside the hot neutron star during the first seconds after the supernova collapse. 
Astronomical observations, such as measurements of the pulsar spindown rate, can be used to 
infer the surface magnetic fields of pulsars some millions of years later. The relation between 
the two is highly non-trivial because of the complex evolution the magnetic field undergoes 
in a cooling neutron star. Let us outline some stages of this evolution. 

Immediately after the formation of the hot neutron star the magnetic field is expected to 
grow due to differential rotation, thermal effects [240], and convection [246]. The dynamo 
effect can probably account for the growth of the magnetic field to about 10^^ — 10^^ G [241]. 

The growth of the magnetic field takes place during the first ten seconds after the supernova 
collapse, in part because the neutrino cooling causes convection. At the same time, during 
the neutrino cooling phase, the neutron star receives a kick. The magnetic field relevant for 
the kick is the average dipole component of the magnetic field inside the neutron star during 
the first 10 seconds of explosion. There is no reason to believe that it has the same direction 
or magnitude as the surface field at the end of the neutrino cooling phase. This, however, is 
only one of several stages in the evolution of the magnetic field. Next, at some temperature 
below 0.5 MeV, the nuclear matter becomes a type-II superconductor. The magnetic field 
lines form the flux tubes, reconnect, and migrate. Next, over some millions of years, the 
pulsar rotation converts the magnetic field energy into radio waves and causes the field to 
evolve even further. The end result of this evolution is, of course, a configuration of magnetic 
fields that is very different from what it was five seconds after the onset of the supernova. 

Clearly, the magnetic field inside a hot young neutron star is not expected to have much 
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correlation with the surface field of a present-day pulsar. 

5.11 Spin-kick from neutrinos 

Spruit and Phinney [204] pointed out that the pulsar rotational velocities may also be ex- 
plained by the kick received by the neutron stars at birth. The core of the progenitor star 
is likely to co-rotate with the whole star until about 10 year before the collapse. This is 
because the core should be tied to the rest of the star by the magnetic field. However, then 
the angular momentum of the core at the time of collapse is 10^ times smaller than the 
angular momentum of a typical pulsar. Spruit and Phinney [204] have pointed out that the 
kick that accelerates the pulsar can also spin it up, unless the kick force is exerted exactly 
head-on. 

The neutrino kick can be strongly off-centered, depending on the configuration of the mag- 
netic field. If the magnetic field of a pulsar is offset from the center, so is the force exerted on 
the pulsar by the anisotropic emission of neutrinos. This mechanism may explain simultane- 
ously the high spatial velocities and the unusually high rotation speeds of nascent neutron 
stars. [200] It was suggested by Phinney (private communication) that a highly off-centered 
magnetic field could be generated by a thin-shell dynamo in a hot neutron star. Since the 
neutron star is cooled from the outside, a convective zone forms near the surface and, at 
later times, extends to the interior. While convection takes place in the spherical shell, the 
dynamo effect can cause a growth in the magnetic field. Thin-shell dynamos are believed to 
be responsible for generating magnetic fields of Uranus and Neptune [247-249]. According 
to the Voyager 2 measurements [250], the magnetic fields of both Uranus and Neptune are 
off-centered dipoles, tilted with respect to the axis of rotation. Unlike other planets, which 
have convection in the deep interior and end up with a well-centered dipole field, Uranus and 
Neptune have thin spherical convective zones near the surface, which explains the peculiarity 
of their dynamos. During the first few seconds after the supernova collapse, convection in 
a neutron star also takes place in a spherical layer near the surface. The thin-shell dynamo 
can, in principle, generate an off-center magnetic field in a neutron star, just like it does in 
Uranus and Neptune. As a result, the anisotropic neutrino emission can give the pulsar the 
kick, and also spin it up. 

5.12 — V correlation 

While the B — v correlation is not expected from a neutrino-driven mechanism, one does 
expect the correlation between the direction of the pulsar motion and the axis of rotation. 
This is because the neutrino anisotropy axis is associated with the dipole magnetic field, 
which is rotating with the star. If the pulsar is accelerated by anisotropic emission of neutrinos 
over the time period of several seconds, the components of the thrust orthogonal to the axis of 
rotation would average to zero, while the component along the axis would not be affected by 
rotation (see Fig. 8). As a result, one expects the pulsar to receive a non- vanishing kick along 
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its rotational axis. To probe this correlation, one needs to perform accurate measurements 
of the polarization of the radio signal from the pulsar. The recent observations confirm the 
fl — V correlation [211,215,216]. In addition, in several cases Chandra images allow one to 
identify the rotational axis of the pulsar with the observed symmetry axis of the pulsar wind 
nebula. These observations also confirm the predicted Q — v correlation [213,218,219]. 




Fig. 8. The anisotropy in outgoing neutrinos, shown here by the density in shading, creates a thrust 
that averages to zero along any horizontal axis, because of the rotation of the protoneutron star 
about the vertical axis. The vertical component of the force, aligned with the axis of rotation, is not 
affected by the rotation. The pulsar velocity should be aligned with the surviving (vertical) compo- 
nent of the neutrino-driven thrust. Hence, one expects a directional correlation between the axis of 
rotation and the direction of pulsar motion. Such a correlation has already been observed [211-219]. 
The rapid motion of the higher-density jet can also be the source of gravitational radiation that 
can be detected in the event of a nearby supernova [16]. 

While the observed correlation between the direction of rotation and the pulsar velocity 
does not provide a definitive prove of the mechanism, it was a generic prediction of the 
neutrino-driven mechanism which turned out to agree with the data. 

5.13 Gravity waves from a neutrino- driven pulsar kick 

In the event of a nearby supernova, the neutrino kick can produce gravity waves that could 
be detected by LIGO and LISA [16,251]. These gravity waves can be produced in several 
ways. 

Obviously, the departure from spherical symmetry is a necessary condition for generating 
the gravity waves. A neutron star being accelerated by neutrinos is not moving fast enough 
to generate gravitational waves from it own motion. However, the anisotropy in the outgoing 
neutrinos turns out to be sufficient to produce an observable signal in the event of a nearby 
supernova. 

Most of the neutrinos come out isotropically and can be neglected. However, a few per cent of 
asymmetrically emitted neutrinos move along the direction of the magnetic field. In general, 
the magnetic field is not aligned with the axis of rotation (Fig. 8), and, therefore, the outgoing 
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neutrinos create a non-isotropic source for the waves of gravity. (Water jet produced by a 
revolving lawn sprinkler is probably a good analogy for the geometry of this source.) The 
signal was calculated by Loveridge.[16] It can be observed by advanced LIGO or LISA if a 
supernova occurs nearby, as shown in Fig. 9. In addition, the neutrino conversion itself may 
cause gravity waves coming out of the core [252-254,251]. 




Fig. 9. Gravity waves signal at LIGO and LISA frequencies calculated by Loveridge [16]. 

The gravitational waves signal from a nearby supernova can be observed by Advanced LIGO 
or LISA. If a fast pulsar emerges from the supernova remnant but the signal is not seen by 
the gravitational waves detectors, it would imply that either the pulsar kick mechanism has 
nothing to do with asymmetric emission of neutrino, or that the magnetic field is aligned 
with the axis of rotation. The radio observations of the nascent pulsar could help determine 
whether the latter is the case. Obviously, the gravitational waves from a nearby supernova 
would be extremely helpful in determining whether the pulsar kick is neutrino-driven or not. 

5.14 Supernova asymmetries 

The supernova explosion can be altered by sterile neutrinos in several ways, in particular, via 
the motion of the neutron star accelerated by the emission of sterile neutrinos. Such a motion 
breaks the spherical symmetry and introduces an additional reason for asphericity that could 
be seen in supernova remnants, as well as optical observations of distant explosions using 
spectropolarimetry and other methods [255-259]. In addition, the overall energy deposition 
behind the shock can increase as a result of convection facilitated by the motion of the 
neutron star [19]. 

Numerical simulations of the neutrino-driven kick mechanism show that the motion of the 
neutron star can seed and drive convection in front of the moving star, adding the energy to 
the ejecta and ultimately helping the explosion [19]. The resultant explosion is asymmetric, 
with the strongest ejecta motion roughly in the direction of the neutron star kick (Fig. 11). 
There should be more mixing in the direction of the neutron star motion, and, in particular, 
nickel should mix further out in this direction [260]. Such an extended mixing in the direction 
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of the neutron star kick has aheady been observed for the fast moving compact star in Sgr 
A East [261]. The asymmetry with an enhancement in the same direction as the neutron 
star motion is in sharp contrast with the consequences of ejecta-driven mechanisms, which 
predict the motion of the ejecta in the opposite direction, in accordance with momentum 
conservation. If both the neutron star and the ejecta receive some additional momentum in 
the same direction, the compensating momentum can be carried by sterile neutrinos. 

The difference between the codirectional asymmetry in the case of a neutrino driven kick 
and counter-directional asymmetry expected from ejecta-driven mechanisms can be used to 
distinguish between the two classes of mechanisms based on the observations of the supernova 
remnants [19]. 

5.15 Sterile neutrinos with the ~ A;eF masses and the energy transport 

Since the sterile neutrinos interact with nuclear matter very weakly, they can be very effi- 
cient at transporting the heat in the cooling proto-neutron star, altering the dynamics of 
the supernova [20,262,263]. In particular, sterile neutrinos with masses in the keV range 
could undergo resonant MSW conversion at some density deep inside the core, below the 
neutrinosphere. This is very different from the active neutrinos, which cannot have an MSW 
resonance at these high densities because their mass differences are too small. The energy 
transfer inside the supernova can change dramatically because of the resonant conversion 
followed by a non-diffusive leap, followed by another conversion, as shown in Fig. 10. This 
could lead to an enhancement of the overall supernova explosion. 
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Fig. 10. The double MSW resonance in the core can facilitate energy transfer out of the core to 
outer regions and can have a dramatic effect on the explosion [20,262] 

An additional enhancement can come from the increase in convection in front of the neutron 
star propelled by the asymmetric emission of sterile neutrinos [19]. Naively one could expect 
that, since sterile neutrinos carry away some energy, the overall energy of supernova should 
be smaller. However, this effect is over-compensated by the enhanced convection in front of 
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the moving neutron star, so that the resulting energy of the shock is higher than in the case 
without sterile neutrinos [19]. 
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Fig. 11. The results of numerical calculations [19] show the increase in convection in front of a 
neutron star moving from left to right due to a neutrino-driven kick. The enhanced convection 
puts more energy behind the shock, which leads to a stronger explosion. Another prediction of 
the neutrino-driven kick is the formation of asymmetric jets, with a stronger jet in the direction 
of he neutron star motion [19]. This is different from hydrodynamic mechanisms, in which one 
expects the stronger jet to point in the opposite direction (due to momentum conservation). Future 
observations may help distinguish between these mechanisms. 



6 Heavy sterile neutrinos and their effect on the supernova explosions 



A substantial range of parameters with large masses and mixing angles (Figs. 1-3) is ruled out 
by the combination of the laboratory experiments and cosmology [79,61]. However, there are 
still windows in which the sterile neutrinos are not only allowed, but also desirable because 
they could augment core collapse supernova shock energies by enhancing energy transport 
from the core to the vicinity of the shock front [264,265,263]. 

One example is a sterile neutrino with mass between 145 and 250 MeV. This range of masses 
is particularly interesting from the point of view of its induced contribution to the mass 
matrix of active neutrinos [6 1] . In this mass range the sterile neutrinos decay predominantly 
into a pion and a light fermion. If the mixing angle with the electron neutrino is negligible, 
and the sterile neutrino mass Ms is in the range m^o < Ms < (m^o -|- m^), the daughter 
pion is the neutral pion, which decays into two photons: [116], where 

a = (/i,T). This decay mode, with lifetime ~ 0.1s, changes the impact of sterile neutrinos 
on the supernova explosion. To distinguish sterile neutrinos that decay mainly into photons 
from the other types, the authors of Ref. [263] called them eosphoric. 

From the ancient Greek god Eujaipdpoi^, the bearer of light. 
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While not in contradiction with supernova 1987A bounds [121], the decays of these neutrinos 
could produce a flux of energetic active neutrinos, detectable by future neutrino observations 
in the event of a galactic supernova [263]. Moreover, the relevant range of sterile neutrino 
masses and mixing angles can be probed in future laboratory experiments [266-270]. 



7 X-ray Detection of Relic Sterile Neutrinos 



The main decay mode of sterile neutrinos in the keV mass range is Vg — > 3i/. This decay mode 
is "invisible" due to the low energy of the daughter neutrinos. In addition to this leading 

mode of decay that occurs through a tree-level diagram, there are also one-loop diagrams 
(Fig. 12) that allow for a photon in the final state. Therefore, the sterile neutrinos can decay 
into the lighter neutrinos and an the X-ray photons: ~^ l^a [271]. The radiative decay 
width is equal to [271,272] 
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and the corresponding lifetime is many orders of magnitude longer than the age of the 
universe. However, since sterile neutrinos are produced in the early universe by neutrino 
oscillations and, possibly, by other mechanisms as well, every dark matter halo should contain 
some fraction of these particles. Given a large number of particles in these astrophysical 
systems, even a small decay width can make them observable via the photons produced in 
the radiative decay. This offers, arguably, the best opportunity to detect these particles. 
Since 1^2^^ — > 71^1™'' is a two-body decay, the resulting photons have energy 

Ey = ms/2, 

which corresponds to a line broadened only by the velocity dispersion of the dark matter 
particles in a given halo. This line, with photon energy of a few keV, can be observed using 
an X-ray telescope [24] . 





Fig. 12. Radiative decay of sterile neutrinos, — > jv^'''' . The X-rays produced by these decays 
can be detected by the X-ray telescopes, such as Chandra, Suzaku, XMM-Newton, and the future 

Constellation-X. 

A broad range of astrophysical systems can provide suitable targets for such observations. A 
concise discussion and comparison of such observational targets can be found in Ref. [102]. 
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First, there should be a signal from distant unresolved sources, in the form of isotropic 
extragalactic background with a (red-shifted) photon line. Second, closer objects with high 
densities of dark matter provide an even better target: the Milky Way halo, the clusters of 
galaxies, and the dwarf spheroid galaxies can give different contributions depending on the 
field of view of one's instrument. In choosing the best target for a specific instrument, one 
usually looks for the most dark matter to fit within the field of view (FOV). This is because 
the contribution of each individual source scales with distance r as while the fraction 

of volume slice of thickness 5r at some fixed distance r that fits within the FOV scales as 
r^. Therefore, to first approximation, it is the total amount of dark matter in the FOV that 
determines the sensitivity of dark matter search. Based on this reasoning, as well as the 
considerations of observational and instrumental backgrounds, dwarf spheroid galaxies are 
probably best suited for searching sterile neutrinos with small-FOV instruments, such as 
Suzaku X-ray telescope [95,96,102]. 

We note in passing that the search strategy for decaying dark matter is quite different from 
the searches carried out for annihilating dark matter, such as super symmetric neutralinos, 
for example. The annihilation rates, and, therefore, the expected signals from annihilating 
dark matter are proportional to the square of the density. Such predictions are sensitive to 
the matter power spectrum on small scales and the morphology of such a signal is difficult 
to disentangle from the point sources. In contrast, sterile neutrino decays are expected to 
produce the fiux proportional to the first power of dark matter number density, not the 
density squared. This fact, combined with the lack of small structure on sub-kpc scales for 
WDM, renders the theoretical predictions for Suzaku [95,96] and other searches fairly robust. 

Several types of constraints have been reported in the literature, and, depending on the 
underlying assumptions, the cxchidcd regions may differ. The X-ray flux depends on the 
sterile neutrino abundance. In view of the variety of possible production mechanisms, one 
can considers two types of X-ray limits that give an answer to one of the following questions: 

• Can a sterile neutrino with a given mass and mixing angle exist (even if it contributes 
only a small fraction to dark matter)? 

• Can a sterile neutrino with a given mass and mixing account for all of dark matter? 

Both questions can be addressed with the use of the X-ray data, but the resulting bounds on 
the mass and the mixing angle can differ significantly. The mass and mixing angle alone do 
not determine the abundance of relic sterile neutrinos. They determine unambiguously the 
amount of relic neutrinos produced by active-to-sterile neutrino oscillations [21], but this need 
not be the only source of sterile neutrinos. Some additional amounts could be produced from 
the Higgs decays [17,33] or the inflaton [65] or radion [150] decays, or any other mechanism, is 
controlled by some additional parameters, unrelated to the mixing angle. It is true, however, 
that, regardless of how many sterile neutrinos were produced at high temperature, the low- 
temperature production by oscillations cannot be switched off (except in the low-reheat 
cosmological scenarios [151]). So, for every mass and mixing angle one can calculate the 
minimal sterile neutrino abundance and determine the number of sterile neutrinos in galactic 
halos using the methods developed in Refs. [44,21-23,25,26,28,30,152,273,63]. The abundance 
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calculated this way gives the lower bound on the sterile neutrino abundance. If this quantity 
exceeds the amount of dark matter determined from cosmological observations, as would 
be the case for large mixing angles, then the particle is ruled out. Otherwise, one can use 
the X-ray observations to set a model-independent bound on the mass and mixing angle for 
the particle that can exist but may contribute only a small fraction of the observed dark 
matter [17,81-83]. Such a particle may not solve the missing matter problem, but it can be 
responsible for the pulsar kicks in some range of parameters. 

Alternatively, one ask whether a sterile neutrino with a given mass and mixing can account 
for all the dark matter, regardless of how it is produced in the early universe. In this case, one 
does not have to calculate the amount produced by neutrino oscillations: the relic abundance 
is assumed to be that of dark matter. This, in turn, determines the flux of X-rays expected 
from a given astrophysical object, for example, a dwarf spheroid galaxy. The X-ray bounds 
obtained this way will, in general, be stronger that those based on production by oscillations 

If all the dark matter is made up of sterile neutrinos {Qg ~ 0.2), then the limit on the mass 
and the mixing angle is given by the dashed line in Fig. 13. However, the interactions in 
the Lagrangian (1) cannot produce such an fi^ = 0.2 population of sterile neutrinos for the 
masses and mixing angles along this dashed line, unless the universe has a relatively large 
lepton asymmetry [22]. If the lepton asymmetry is small, the interactions included in the 
Lagrangian of eq. (1) can produce the relic sterile neutrinos only via the neutrino oscillations 
off-resonance at some sub-GeV temperature [21]. The model-independent bound [17,81-83] 
based on this scenario is shown as a solid (purple) region in Fig. 13. It is based on the flux limit 
from X-ray observations [24,97-110,103,111,112,122] and the state-of-the-art calculation of 
the sterile neutrino production by oscillations [30,152]. 

Finally, there are differences in the way the data are analyzed in search of the line from the 
sterile neutrino decay. For example, some limits are based on the requirement that the X-ray 
flux from dark matter not exceed the observed total flux from a given object. Some other 
limits are derived using the background modeling based on known astrophysical sources. 
Finally, yet another strategy is to use an ad hoc numerical fit to the data and to require that 
the goodness of fit not be affected by the addition of a line at a given energy. While each of 
these strategies can be defended, it is obvious that the former approach produces the most 
conservative limits, while the latter produces the strongest, but less robust limits. 

Instrumental backgrounds and energy resolution often limit one's ability to detect a line. 
While some existing telescopes, such as, e.g., Suzaku have a relatively stable background 
and reasonable energy resolution, the advent of new instruments, such as Astro-H and 
Constellation-X, with a much better energy resolution will allow a much more efficient search. 
Meanwhile, the ongoing search using Suzaku X-ray telescope targets the dark matter dom- 
inated dwarf spheroid galaxies, such as Draco and Ursa Minor [95,96]. The outcome of this 
search may allow one to explore a substantial part of the range of parameters in which the 
pulsar kicks and dark matter could be explained simultaneously. 
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Fig. 13. The solid excluded region is based on a combination of the X-ray and cosmological bounds; it 
applies even if sterile neutrinos constitute only a fraction of dark matter [17,81,122]. The dashed line 
shows the X-ray bound under the assumption that sterile neutrinos (produced by some mechanism 
that may be different from oscillations) make up all the dark matter. Additional bounds from 
small-scale structure formation (see discussion of Lyman-a observations) may apply, depending 
on the free-streaming length, whose relation with the particle mass depends on the production 
scenario [157,156,167,82,83,184]. The allowed regions for the pulsar kicks shown here are based on 
Refs. [11,18,13]. 

8 Early decays, X-rays, and the formation of the first stars 



The X-ray photons from sterile neutrino decays in the early universe could have affected the 
star formation. Although these X-rays alone are not sufficient to reionize the universe, they 
can catalyze the production of molecular hydrogen and speed up the star formation [51- 
54], which, in turn, would cause the reionization. Molecular hydrogen is a very important 
cooling agent, necessary for the collapse of primordial gas clouds that gave birth to the first 
stars. The fraction of molecular hydrogen must exceed a certain minimal value for the star 
formation to begin [274]. The reaction 



H + H^H2 + 7 



(52) 



is very slow in comparison with the combination of reactions: 



H 
H 



Ho + H^ 



(53) 
(54) 
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which become possible if some hydrogen is ionized. Thns. the ionization fraction determines 
the rate of molecular hydrogen production. If dark matter is made up of sterile neutrinos, 
their decays can produce a sufficient flux of photons to cause a signiflcant increase in the 
ionization fraction [51 54]. This can have a dramatic effect on the production of molecular 
hydrogen and the subsequent star formation. 

Decays of the relic sterile neutrinos during the dark ages could produce an observable signa- 
ture in the 21-cm background [275]. It can be detected and studied with such instruments 
as the Low Frequency Array (LOFAR), the 21 Centimeter Array (21CMA), the Mileura 
Wide-field Array (MWA) and the Square Kilometer Array (SKA). 



9 If this is right... 

If, indeed, the relic sterile neutrinos are discovered, the sharp spectral line from their two- 
body radiative decay will be a boon to observational cosmology, because one can hope to 
get the distance-redshift information from the X-ray observations of dark matter halos. In 
addition, the X-ray instruments with high energy resolution and good imaging capabilities 
could be used to map out the dark matter halos that contain no luminous matter. 



10 Conclusions 

The fundamental physics responsible for the neutrino masses is likely to involve some gauge- 
singlet fermions. Their Majorana masses can range from a few eV to some values well above 
the electroweak scale. If some of the Majorana masses are below the clectroweak scale, the 
corresponding new degrees of freedom can have direct observable consequences. In particular, 
a sterile neutrino with mass of several keV is a viable dark matter candidate [21], and the 
same particle can explain the pulsar kicks [13,11] and can play a role in other astrophysical 
phenomena. 

The search for an X-ray line from the radiative decay of relic sterile neutrinos provides the 
best opportunity to discover these particles if they exist. Depending on the model parameters, 
relic sterile neutrinos can constitute all or part of cosmological dark matter. Therefore, dark 
matter dominated systems, such as dwarf spheroid galaxies, which are also X-ray quiet, are 
among the best targets for current and planned X-ray telescopes. In addition, the diffuse 
extragalactic background and the Milky Way halo can be used to search for sterile dark 
matter. The improved energy resolution of the next generation of X-ray telescopes, such 
as Astro H (know as NeXT in its earlier stages of development), will help explore the best 
motivated region of parameters, in which both dark matter and the pulsar kicks are explained 
by the same sterile neutrino. Due to the small mixing angles, the laboratory experiments 
will probably not be competitive with the astrophysical searches in the foreseeable future. 
However, if sterile neutrinos are to be a relatively cold dark matter, the production of these 
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particles should have taken place at or above the electroweak scale, and LHC may be able 
to discover the requisite gauge singlet in the Higgs sector [17,33]. The observations and the 
population studies of pulsars can provide some additional clues regarding the role of sterile 
neutrinos in generating the pulsar kicks. Given the broad range of X-ray searches that have 
been conducted and are being planned, one can hope to explore the best-motivated range of 
masses and mixing angles in the near future. 
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